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Fluid and solute transfer in X-junctions between straight channels is shown to depend critically on 
the junction angle a in the Stokes flow regime. Experimentally, water and a water-dye solution are 
injected at equal flow rates in two facing channels of the junction: Planar Laser Induced Fluorescence 
(PLIF) measurements show that the largest part of each injected fluid "bounces back" preferentially 
into the outlet channel at the lowest angle to the injection; this is opposite to the inertial case and 
requires a high curvature of the corresponding streamlines. The proportion of this fluid in the 
other channel decreases from 50% at a = 90° to zero at a threshold angle. These counterintuitive 
features reflect the minimization of energy dissipation for Stokes flows. Finite elements numerical 
Ch . simulations of a 2D Stokes flow of equivalent geometry confirm these results and show that, below 

the threshold angle ac = 33.8°, recirculation cells are present in the center part of the junction and 
separate the two injected flows of the two solutions. Reducing further a leads to the appearance of 
new recirculation cells with lower flow velocities. 



Flow control in mierochannels has become an important area of research in microprocess engineering^. 
The behavior of the fluids at junctions in microfluidic circuits is particularly critical for applications such 
■ as mixing, chemical reactions or heat exchange: this is for instance the case when two different fluids are 

injected separately into a junction of flow channels. While, for T-junctions, the flow of the two fluids at the 
outlet depends weakly on the junction angle at low Reynolds numbera^i^, there are very diverse patterns in 
X-junctions. Then, the distribution of the injected fluids between the two outlets and their interaction within 
the junction, which are both crucial for applications, depend strongly on the junction angle a (0 < a < 90°). 

For this reason, while most previous studies dealt only with orthogonal^ or parallel^i^ channels, the 
present work analyzes specifically the infiuence of the angle a. More precisely, the two fluids are injected 
in facing channels and we study as a function of a their relative fraction in each of the outlet channels 
and the appearance of viscous eddies. Other authors used tangent channels at a varying angled but the 
interaction between the streams differs strongly from the present case of channels in a same plane. This 
system is studied by combining experimental PLIF measurements on transparent models and 2D numerical 
simulations. Like for Fan and Hassan^, the experimental system is upscaled in order to visualize conveniently 
the flow; however, the Reynolds number Re = paU/rj is low enough so that inertial effects are negligible and 
the Stokes equation is satisfied (U, p and rj are respectively the mean velocity, density, and viscosity of the 
fiuid and a is the width of the channels). The flow velocity field is therefore the same as in much smaller 
mierochannels of same geometry. 

The channels are carved into a transparent plexiglas plate by a computer controlled milling machine; 
both their width a and depth h are constant and equal to 3.55 mm. The horizontality of the plate is 
carefully adjusted. Each outlet is connected to one of the inlets of a double syringe pump sucking the 
fluids at equal flow rates Qoi and Qo2 (FigH]). The inlet channels iw and id are connected to two glass 
flasks containing respectively pure water and a solution of water and rhodamin 6G dye (0.25 mg.l~^). Due 
to this low concentration, the two fluids are Newtonian: their viscosities and densities are so close when 
their temperatures are matched that no buoyancy or viscosity driven instabilities may appear. Each bottle 
is placed on computer controlled scales allowing for a measurement of the flow rates Qi^ and Qid in the 
corresponding outlet. In the present work, these flow rates are set to be equal within ±0.3% by adjusting 
the relative levels of the fluid in the flasks. 
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FIG. 1. Experimental setup. 



The distribution of the two fluids in the intersection is visuaUzed by the Planar Laser Induced Fluorescence 
(PLIF) technique. A plane horizontal green laser sheet of wavelength 532 nm and thickness ~ 1 mm, parallel 
to both channels, enters the plexiglas model from the side. The parts of the channels containing the rhodamin 
dye solution are easily identified by the emission of a yellow fluorescence light of intensity proportional to 
the local dye concentration. The distribution of this solution in the illuminated plane is then recorded 
by a digital camera located above the model : unwanted reflected light from the laser is eliminated by a 
narrow bandwidth notch filter corresponding to its wave length A = 532 ± 25 nm and placed in front of the 
camera sensor. Regions in which only water is present appear as dark while those containing rhodamin are 
illuminated (Fig. [2]) . In the present work, our analysis is based on the geometry of the boundary between 
the zones occupied by the two fluids. The laser is mounted on a vertical translation stage: this allowed us 
to scan the distribution of the fluids at different heights z for several Reynolds numbers (0.5 < Re < 50). 

The distributions of the fluids in the symmetry plane at different Re values have first been compared 
in order to check the validity of the linear Stokes equation. For Re — 0.5 and Re — 5, the geometry of 
the boundary between the transparent and dyed fiuid is the same: this implies that the linearity condition 
is satisfied and that the fiow is purely viscous. At higher Re values {Re > 15), inertial effects induce, 
as expected, changes of the shape of the boundary. Moreover, for Re = 0.5 and Re = 5, there is no 
visible influence of the vertical coordinate z of the light sheet: the distribution of the two fluids can then be 
considered as bidimensional. The flow fleld itself is however not bidimensional because of the vertical velocity 
gradients created by the upper and lower wall: however, the orientation of the velocity and, therefore, the 
streamlines are approximately invariant with z. In contrast, for Re > 40, 3D structures appear in the fluid 
distribution. The boundary between the dyed and the transparent fluid in the images is fuzzier for Re = 0.5 
than for Re = 5 and the thickness of the transition zone increases with distance to a value of the order 
of 0.5 mm for a path length of the order of a cm . This value is of the same order of magnitude as that 
corresponding to transverse molecular diffusion (using Dm — 2.8 lO^^'^ m^.s~^ at T = 22°C^). 

Taking into account the above results, all experiments reported below were performed at Re = 5 for 
which both sharp boundaries and Stokes flow are obtained. Cells with 6 different junctions angles a = 
90,80,45,30,20 and 10° have been used. Figure [5] displays pictures obtained in the stationary regime for 
three angles a = 80°, 45° and 30°. The dark and illuminated regions of the junction are symmetrical with 
respect to the center of symmetry of the junction : this confirms that the fiow rates Qid and at the 
two inlets are equal. A first key observation is that, for all angles a, each injected fluid flows dominantly 
towards the outlet at the lowest angle with the injection channel (outlet 1 for the fluorescent solution): this 
corresponds to flow paths with a higher curvature than those leading to the other outlet. The flow of each 
fluid is equally split between the outlets only for a = 90°. At a = 45°, the fraction of dyed fluid moving 
towards outlet 2 is small and is exactly zero for a = 30°. In this latter case, all incoming streamlines "bounce 
back" on the junction towards outlet 1. The main flows of the two fluids are then completely separated. 

More quantitatively, the distribution of each fluid between the two outlets will be characterized by the 
volume fraction q occupied by the fluid of interest: practically q is taken to be equal to Ud/a m which a is 
the total width of the section and ad the width of the flow tube corresponding to this fluid in the parallel 
flow region. Here, ad is computed for the dyed solution in outlet 2 where its volume fraction is lower (Figs. [T] 
and[2^). Several measurements of ad in the parallel flow region are averaged in order to compute q. Figure|3] 
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FIG. 2. Experimental images displaying the distribution of the fluids in the junction for Re = 5. (a) : a = 80°), (b): 
a = 45°, (c): a = 30°. 




FIG. 3. Variation of ratio q = ad /a (see Fig. [2^) as a function of the junction angle a. (•) : experimental 
measurements; (+) : results of the numerical simulations. 



displays the experimental variation (• symbols) of the fraction q with a: q is exactly zero below a threshold 
angle ac of experimental value in the range 30° < olc < 35°. Above etc, 9 increases linearly with a up to 
q = 0.5 for a = 90°. 

The zero value of q below the critical angle ac may seem at first counter-intuitive: it implies indeed that the 
fluid particles follow preferentially, in the region of the junction, the streamlines with the highest curvature. 
Actually, this reflects a key characteristic of all viscous Stokes flows, namely that they must achieve a 
minimum viscous energy dissipation. Let us assume that, in contrast to the experiment of Fig. [5};, the dyed 
fluid flows into outlet 1 and the transparent one into 2 at low values of a. This requires a counterflow of the 
two fluids along segment CD which is much longer than AB] moreover, this small length of AB increases 
the transverse velocity gradient. The viscous energy dissipation, proportional to the integral of the square 
of this gradient, would then be too high. In the experimental flow of Fig. [5J:, instead, the distance AB along 
which the two flows coexist is smaller than the transverse distance CD: together with the appearance of 
recirculation cells with very low velocities (and gradients), this leads, indeed, to a lower dissipation than in 
the first case. The result would be completely different for the inertial potential flow of a perfect fiuid in 
the same geometry. Let us consider, for instance, the comparable classical problem of a 2D jet of perfect 
fluid impinging obliquely at an angle a on a plane walli^i: then, as a decreases from 90° to zero, the fraction 
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(equivalent to q) of the fluid following the path with the lowest curvature increases from 0.5 to 1 (instead 
of decreasing to like in Fig. [3]). The structure of the flow in the center part of the junction is obviously a 
key element for predicting the exchanges between the two flowing fluids. Additional information is obtained 
by saturating the intersection with the rhodamin- water solution instead of pure water before establishing a 
stationary flow in the facing injection channels. No difference is observed for a > Uc but, for a < ac, the 
fluorescence light intensity is higher in the center of the junction than when it was initially saturated by 
transparent fluid. This suggests that some dye is trapped in the center of the junction and that the latter 
corresponds to a dead or recirculation zone. 




FIG. 4. Velocity field (vectors) obtained numerically for : (a) a = 40°; (b) a = Oc = 33.8° (c) a = 31°. The length 
of all vectors are equal for a better visibility so that only the orientation of the velocity is represented. Continuous 
lines represent selected streamlines of the flow; dashed lines mark the boundaries of the flow recirculation regions. 



In order to verify this assumption, the 2D Stokes equation has been solved numerically in the (x, y) plane 
by means of the FreeFem finite elements package^^ : the lack of inertial effects and the 2D distribution of 
the two fluids suggest indeed that this equation may predict adequately the distribution of the two fluids. 
The velocity vector field and the sample streamlines obtained in this way for three values of the angle a 
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FIG. 5. Typical streamlines obtained numerically for different angles (a): q = 35°), (b): a = 24°, (c): a = 19°, (d): 
a = 15°. Continuous and dashed lines have the same meaning as in Fig. |4] 

are displayed in Fig. 2] For a = 31° (Fig. |4t), a recirculation cell is observed in the center of the junction 
as suggested by the experiments. It separates completely the flows on the two sides of the junction. The 
boundary of this zone inside which the flow lines are not connected to the outside is marked by dashed lines. 
The modulus of the velocity (not represented on the figure) is about 100 times lower inside the cell than 
outside. The detailed structure of the recirculation flow depends on the angle a and will be discussed below. 
For a = 40° (Fig. SK), there is only one stagnation point and it is located at the center of symmetry of the 
junction: streamlines originating at this point separate those extending across the junction from those which 
"bounce back" on it. The distance of these particular streamlines to the nearest wall represents therefore 
the distance ad defined above. This provides a numerical prediction of the ratio q = ad/ a: its variation as a 
function of the angle a is displayed by (+) symbols in Fig. [H The numerical values of q are first very close 
to the experimental ones: this supports (a posteriori) the 2D approximation. The variation of q with a is 
observed to be very linear except close to the critical angle ac- this confirms the trend of the experimental 
data. These simulations provide the precise value ac = 33.8 ± 0.2°. For a = ac (Fig. Hb), one observes 
two recirculation cells on the lateral edges extending up to the central stagnation point. For this and lower 
angles, there is no fiow path connecting the left and right sides of the junction any longer. 

The recirculation zones observed at low angles a clearly play a major part in the exchange of passive scalars 
(solutes, heat . . .) between the injected flows. Their structure depends strongly on a as shown in Figure [5] 
displaying typical examples obtained numerically. For a slightly above ac (Fig. [5^), side recirculation cells 
are present, like for a = ac'. they reduce the section available in the plane of symmetry (vertical on the 
flgure) for streamlines crossing from one side to the other but do not block it. This section becomes zero for 
a = ac as the two cells touch each other (Fig.Hb) and then, below ac, a single cefl extending all across the 
junction appears (Fig. [41:). The faster variation of q near ac in Fig. |3]may reflect the scarcity of the flow 
paths connecting the two sides of the junction and the fact that they cross low velocity regions. Decreasing 
further a (Fig. \E}p) results in an elongation of the cell in the x direction with, inside it, two smaller cells 
joined at a stagnation point (Fig. [5b)- At a still lower angle, the latter occupy nearly the full width of the 
channel (Fig. [St) until an additional cell, elongated in the y direction, appears in the center (Fig. [SJi). The 
typical velocity decreases by a large factor (~ 100) between the outer flow and the first generation of cells 
and between the latter and the next generation. New cells appear at still lower angles (not shown on the 
figure) and correspond to vanishingly small velocities. 

Vortices have already been observed in Stokes fiow by many other authors. A classical example is Moffatt 
vortices induced by an outside flow in a corner of angle f3 < 146° between plane walls: an inflnite sequence 
of eddies appears with velocities decreasing strongly, like here, from one generation to the nex t^^i^^ . Here, a 
similar sequence appears at low angles a but is interrupted at the middle of the junction with a symmetrical 
sequence starting on the other side. These Stokes flows in corners do not display however like in Figs. [S^-c 
pairs of corotating vortices separated by a stagnation point unless a planar Couette flow is superimposed^^. 
These latter structures are also observed in elongated rectangular cavities when the two smaller parallel 
walls are sheared in opposite directional^. Like at low a values in the present system, as the elongation of 
the cavity increases, new recirculation cells appear : one observes flrst a pair of vortices and they coalesce 
into a single cell when the elongation is increased further. 

The experiments and simulations presented here demonstrate that, at all a values (except 90°), the largest 
part of (if not all) each injected flow bounces back into the outlet at the angle a to the injection. The actual 
fraction depends strongly on a, as shown above, but will also vary with the aspect ratio h/a oi the section: 
for instance, if /i/a -C 1 (instead oi h/a = 1), one obtains a Hele Shaw cell geometry. Then, the average 
of the velocity field over the depth h would be similar to that for a perfect fluid and the variation of q 
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with a will be completely different. It will be important to study the transition between these two opposite 
behaviors by varying continuously h/a. Since the flow field is controlled by the viscous dissipation, it is 
likely to be strongly influenced by the rheological properties of the fluid. Moreover, anomalies are known 
to occur for fluids of high cxtcnsional viscosity in the vicinity of the stagnation points in different flow 
geometriesi^"— . For viscoelastic fluids flowing in orthogonal channels, elastic instabilities have indeed been 
observed experimentally^^ and flow asymmetries have been predicted numerically^!. 

The recirculation cells appearing at low a values will control the exchange of solute (or other passive or 
reactive species) between the two injected fluids: this exchange involves a combination of transverse molec- 
ular diffusion between the cells and the two flows (and from one cell to another) and convective transport 
by the circulation within the cells. For reactive fluids, the recirculation cells may act as microreactors into 
which species contained in the injected fluids are transferred through transverse molecular diffusion. It may 
be possible to induce a more efficient chaotic mixing in the recirculation cell by varying the geometry of the 
intersection and/or by using suitable time variations of the relative values of the two injection flow rates. 
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